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A vortex in a model spinless px -\-ipy superconductor induces two Majorana fermions (MFs), one 
in the core and the other at the sample edge. In the present work, we show that edge MP can be 
generated, fused, transported, and braided easily by tuning gate voltages at point-like constriction 
junctions. Solving the time-dependent Bogoliubov-de Gennes equation, we demonstrate that the 
braiding of edge MFs obeys the non-Abelian statistics. The present setup is therefore a promising 
implementation for topological quantum computation, and has the advantage of easy manipulation 
and simple device structure. 

PACS numbers: 



Topological quantum computation is attracting consid- 
erable interests due to the unique feature of fault toler- 
ance, where quantum information is stored non-locally 
and robust against decoherence caused by interaction 
with environment [ll, [2| . Ground state degeneracy and 
rotations within the degenerate subspace with unitary 
non-Abelian transformation are the two important ingre- 
dients for its implementation. Systems with Majorana 
fermions (MFs) |3l42^, particles equivalent to their an- 
tiparticles, as zero-energy excitations are promising, since 
they are half of conventional fermions and can form non- 
local qubits. Superconducting state provides a natural 
host for these zero-energy MFs where Bogoliubov quasi- 
particles are composed by both electrons and holes. 

It was illustrated that a vortex in spinless Px + Wy 
superconductor can accommodate a zero-energy MF at 
its core [5[, and the braiding of vortices obeys the non- 
Abelian statistics [6]. A hetero structure of 5- wave su- 
perconductor (S) and topological insulator (TI) (25l-[27j 
was then proposed to behave similarly to the spinless p- 
wave superconductor due to strong spin-orbit coupling, 
and circuits of S/TI/S junctions can achieve creation, 
manipulation and fusion of MFs by tuning the super- 
conductivity phases [11]. Very recently, it was demon- 
strated that a spin-orbit coupling semiconductor (SM) 
in proximity to a ferromagnetic insulator (FI) plays a 
similar role as TI, and thus a S/SM/FI hetero structure 
can be a generic platform to provide topological phase 
It was also illustrated clearly by a toy model that 
a one-dimensional (ID) spinless superconductor carries 
two zero-energy MFs at its two ends [7], which can be 
realized by a ID spin-orbit coupling semiconductor un- 
der a magnetic field and in proximity to an 5-wave su- 
perconductor [isl, [l9| . The non-Abelian braiding of end 
MFs has been demonstrated very recently in a network 
of these nanowires with voltage applications (2Q| . 

In the present work, we notice that the edge MF real- 
ized in a finite sample with a superconducting vortex is 
very useful. Since MFs always come as pairs, the zero- 
energy MF bounded at the vortex core has its counter- 
part, which appears at the sample edge (see Fig. 1(a)). 




FIG. 1: (a) Schematic device setup for generating and braid- 
ing of edge MFs by tuning gate voltage at point-like constric- 
tion junctions, (b) Each of the four finite samples (called 
bricks) consists of a hetero structure of s-wave superconduc- 
tor with a vortex, spin-orbit coupling semiconductor and a 
ferromagnetic insulator, (c) Tuning chemical potential to the 
gap between the two bands, a superconducting gap is opened 
at the outside one. 



The edge MFs have important advantages in manipu- 
lation, namely they can be created, fused, transported 
and braided easily by applying gate voltages on constric- 
tion junctions between finite samples as revealed in the 
present work. In contrast, motion of core MFs should 
always be accompanied by vortices, which is not an easy 
task practically. Comparing with previous proposals 
based on control of superconductivity phases in circuit 
of S/TI/S junctions [11| and manipulation of end MFs in 
network of ID nanowires [20], gate voltages are applied at 
the point-like constriction junctions in the present setup, 
which makes the device structure and operations sim- 
pler. As a drawback of the present device, the energy 
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gap between edge zero-energy MFs and finite-energy ex- 
citations is smaller compared with the core counterparts, 
thus requiring manipulation at lower temperature. It is 
expected, however, that this difficulty can be overcome 
relatively easier by using superconductor with large en- 
ergy gap. 

Explicitly we consider a finite sample of the S/SM/FI 
hetero structure with a superconducting vortex at the 
sample center 0, which we call a brick. Bricks are 
connected by constriction junctions, where gate voltages 
can be tuned to connect and disconnect the constriction 
junctions adiabatically. We show that edge MFs can be 
created, fused, transported, and braided with sequences 
of switching on and off the gate voltages in the sys- 
tem. Solving the time-dependent Bogoliubov-de Gennes 
(BdG) equation, we monitor the time evolution of the 
wave functions, and demonstrate clearly that braiding of 
edge MFs obeys the non-Abelian statistics. 

Bogoliubov-de Gennes Hamiltonian 

Our elementary system is schematically depicted in 
Fig. 1(b). The Hamiltonian of a spin-orbit coupling semi- 
conductor in proximity to a ferromagnetic insulator is 
given by 



Ho = J drc\r) 



2m* 



c(r) 



(1) 

with m*, /i, aR and Vz being the effective electron 
mass, chemical potential, strength of the Rashba spin- 
orbit coupHng, and Zeeman field, respectively, and a = 
{cFx^cFy^cFz) the Pauli matrices and c = (c^ c^)^ the elec- 
tron annihilation operators. The proximity effect from 
the 5-wave superconductor is described by 



where A(r) is the effective s-wave pairing potential. The 
total Hamitonian is i^tot = + ^sc- For chemical po- 
tential Y^A^ + /i^ < Vz , the system behaves effectively as 
a 2D spinless Px^Wy superconductor with the dispersion 
and spin configuration given in Fig. 1(c) |5l, fiol, l28l - [3Qj . 

The quasiparticle excitations in the system are de- 
scribed by the following Bogoliubov-de- Gennes (BdG) 
equation 



A 

At -CFyH^CFy 



^(r) =^^(r). 



(3) 



where the Nambu spinor notation ^(r) = 
[u^{f)^ui{f)^vi{f)^—v^{r)\ is adopted. Since Hq is 
Hermitian, any eigen vector ^ of Eq.(3) with energy E 
has its twin cFyiy^* with energy —E^ where Ty is the 
Pauli matrix in Nambu spinor space. Because the BdG 
equation is defined by an even dimensional Hamiltonian, 
zero-energy eigen modes, if any, should appear in pairs. 
By recombining these zero-energy eigen wave functions. 



one can always achieve ^ = a^yf^y^*, for which the 
quasiparticle operator defined by 



(4) 



satisfies the relation 7' = 7. Therefore, the zero-energy 
excitations of the system are actually pairs of MFs. 

It has been shown [10, 29, 30] that the system Eq. (3) 
with a vortex in the superconductor accommodates a MF 
at the vortex core. For the case that there is only one 
vortex in the s yst em, its twin should appear at the edge 
of the system [2^, [3l| . While the core MF has been in- 
vestigated in many previous studies, little attention has 
been paid on the edge MF (see [11]), which is the main 
focus of the present work. 

In order to explore the edge MF, we resort to numerical 
analysis. For this purpose, we derive the tight-binding 
version of Hq on a square grid. 
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where both spin-reserved hopping to = /2m*o? and 
spin- flipped hopping ^q, = aRh/2a are between nearest 
neighbors with a the grid spacing, and ji is measured 
from band bottom. 

The tight-binding BdG equation has the same form of 
Eq.(3), and becomes AN coupled linear equations with 
N being the total number of grid sites. The energy spec- 
trum and all wave functions of excitations including those 
of MFs can be obtained by diagonalizing the 47V x AN 
matrix. In the present study, the square sample is di- 
vided into a grid of 100 x 100 sites. Although the size 
of the Hamiltonian matrix can be large typically in the 
order of 10^, in most cases only several states near E = {) 
are relevant^ for which a powerful numerical technique is 
available [3/ 

Edge Majorana Fermions 

We first study a square sample (see Fig. 1(b)) with a vor- 
tex at the center A(r) = Ao[l - e-(^/^)']e^^ with s = 4a, 
namely the brick. As shown in the energy spectrum in 
Fig. 2(a), we found a pair of zero-energy excitations in 
the system. By recombination, we get two wave func- 
tions, one at the vortex core and the other on the sample 
edge with their norms plotted in the inset of Fig. 2(a), 
both satisfying the generic relation for MF. We define 
them as core MF and edge MF. This treatment of wave 
functions is understood for all discussions through the 
present work. 

Next we study two bricks connected by a constriction 
junction. As shown in Fig. 2(b), we find again two zero- 
energy MFs, each bounded by one of the two vortices. 




1 2 3 4 5 6 

n 



FIG. 2: Energy spectrum of several lowest excitations and 
distribution of wave-function norms of zero-energy MFs, (a) 
for one brick with two MFs, one at vortex core and the other 
at brick edge, and (b) for two connected bricks with two core 
MFs. Results are for Ao = 0.5to, Vz = 0.8to, M = and 
to, = 0.9to, with 100 X 100 sites for the brick and 10 x 10 sites 
for the constriction junction. 



There is no edge MF this time (see inset of Fig. 2(b)), 
since the two edge MFs in the two bricks meet each other 
at the constriction junction and fuse into a Bogoliubov 
quasiparticle with finite energy. 

We proceed to investigate a system of three bricks, 
and introduce dynamic process of switching on and off 
constriction junctions. This is implemented by tuning 
the hopping integrals at the constriction junctions among 
bricks, which can be realized practically by adjusting gate 
voltage. At the initial state, the left brick is isolated 
while the other two are connected. There are three core 
and one edge MFs as revealed above. As confirmed by 
calculation, the core MFs remain unchanged (in an expo- 
nential accuracy) during the switching process, we omit 
them in the following discussions for simplicity. The con- 
striction junction between the left and central bricks is 
then turned on adiabatically. As displayed in Fig. 3(a), 
it is interesting to observe that, the edge MF on the left 
brick spreads its wave function to the central and right 
bricks, since now there is only one unified edge of the 
three bricks. We emphasize that this is one of the clear- 
est manifestations of the mobility of the edge state. 

We then turn off the constriction junction between the 
central brick and the right one adiabatically. As shown 
in Fig. 3(a), the wave function of edge MF now shrinks 
to the right brick. As the result of these two switchings, 
the edge MF initially at the left brick is transported to 
the right brick. 

During the whole process, we monitor the energy spec- 
trum and make sure that the gap between the zero-energy 
MFs and excitations of finite energy remains open (see 
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FIG. 3: (a) Diagram for transport of edge MF prepared at 
the left brick to the right brick and distributions of calcu- 
lated wave-function norms, (b) Evolution of the energy spec- 
trum, where the two zero-energy excitations are associated 
with three core and one edge MFs, when the left junction 
is turned on adiabatically parameterized by Ai, and then the 
right junction is turned off adiabatically parameterized by A2. 



Fig. 3(b)), which protects the topological phase of the 
system. Compared with previous prescriptions jll|, [20| , 
the present manipulations on edge MFs are performed by 
applications of point-like gate voltages, which makes the 
device simpler. 

Braiding with Non-Abelian Statistics 

Based on the above results, we can exchange edge MFs 
pairwisely. The simplest structure for this purpose con- 
sists of four bricks, as depicted in Fig. 1(a), and the di- 
agram of exchange process of two edge MFs are demon- 
strated in Fig. 4. As the initial state (step-0), there are 
two edge MFs located at the left and right bricks, whereas 
the central and top bricks are connected (thus no edge 
MFs). Following the process for transportation revealed 
above, we first transport the red edge MF to the top brick 
(step-1, -2). We then transport the green edge MF from 
the right brick to left brick (step-3, -4). Finally, we bring 
the red edge MF tentatively stored at the top brick to 
the right brick (step-5, -6). After these operations, the 
setup of the system, or the Hamiltonian, returns to the 
initial one, but leaving the two edge MFs exchanged. 

Now let us reveal the impact of this exchange to the 
wave function of system. For this purpose, we evalu- 
ate the time evolution of wave function upon adiabatic 
switchin gs by using the time-dependent BdG (TDBdG) 
method [33-[35|. We first derive the wave functions of 
edge MFs at the left and right bricks by diagonalizing the 
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FIG. 4: Diagram for exchanging two edge MFs (red and green for eye-guide) and distribution of wave-function norms. 
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FIG. 5: Time evolution of the wave function of edge MFs for 
the exchange process in Fig. 4, in terms of its projections Ol 
and Or to the two edge MFs prepared at the initial stage. 
The time for a transport process is taken as T = 40000/i/to. 

Hamitonian Eq.(3), with which the initial wave function 
for the total system is given by = 0)) = |(/>l) + \ 4'r)' 
Upon on the switching process illustrated in Fig. 4, the 
wave function evolves with time according to the TDBdG 
equation: 



in^^ = H^ot{m{t)). (6) 

The evolution of wave function can be monitored by its 
projections to the initial wave functions of the two edge 
MFs, Ol = andOi? = At the initial 

state (step-0), one has Ol = Or = +1 by definition. As 
shown in Fig. 5, the order parameter Ol drops gradually 
since the wave function of the red edge MF spreads to the 
three bricks when the left junction is turned on (step- 
1 in Fig. 4); it drops to zero when the red edge MF 
is totally transported to the top brick (step-2). When 
the right junction is then turned on and thus the green 
edge MF spreads its wave function to the central and left 
bricks (step-3), Ol becomes finite, and most intriguingly 
it acquires a minus sign. It achieves Ol = —1 when 
the left junction is switched off (step-4), and remains its 
value for the final two steps. 

The evolution of the order parameter Or up to step-4 
is straightforward. When the red edge MF stored tenta- 
tively at the top brick spreads to the right brick (step-5 



and -6), Or becomes positive and then achieves Or = +1 
by the time the whole process is finished. 

Therefore, during the process of exchanging the two 
edge MFs illustrated in Fig. 4, the wave function of the 
system evolves as 

= 0)) = \4>l) + \<Pr) ^ = 6r)) = -\<Pl) + \4>r), 

(7) 

or in terms of operators, 

7l IR, in -7l, (8) 

which can be presented by a unitary transformation 
Ulr = 6 4^^^^, the same form observed first for core 
MFs [6] (see also [20^). 

The above exchange rule of edge MFs can be under- 
stood in the following way. Let us start from the MF 
transport in Fig. 3. An effective Hamiltonian for the 
low-energy physics of the adiabatic transport is given by 

^eff(t) = ^Ai(t)rLC7L7C + i>^2it)^CRlClR, (9) 

where Vij = —Tji denotes interaction between MFs, and 
for edge MF transport \i(t) increases from to 1 from 
t = to t = T while A2 remains unity, and from t = T 
to t = 2T X2{t) decreases from 1 to while Ai remains 
unity. A zero-energy edge MF can be composed from the 
three MF creation operators at any moment provided the 
topological phase is protected by the adiabatic processes 

iMt) , lR\l{t)TLc/TcR 

(10) 

with 7(t = 0) = 7l and ^{t = 2T) = sgn{TLC /^CR)7R' 
An interaction between jl and jr can be included, but 
it is easy to see that the result of transport remains the 
same. 

The exchanging process which consists of three trans- 
port processes with two edge MFs as shown in Fig. 4 
is then described in the following way: 71 (t = 2T) = 
sgn(rLc/rcT)7T and 72 (t = 2T) = 7^^, 71 (t = 4T) = 



5 




4 6 8 

time [ T ] 



10 12 



O 



0.0 




10 12 



time [ T ] 



FIG. 6: Diagrams for two braiding processes of three edge 
MFs (red, green and blue for eye-guide) and time evolution 
of the wave function of edge MFs, in terms of its projections 
Ol, Or and Ot to the three edge MFs prepared at the initial 
stage. 



7i(t = 2T) and j2{t = 4T) = sgn(r^c/rcL)7L, 
71 = 6T) = sgn{TTc/rcR)^gn{TLc/rcT)lR = 
-sgn{r Lc/rcR)lR and j2{t = 6T) = 72(t = 4T). It 
is then clear that after the braiding, the two MFs at the 
left and right bricks acquire opposite signs, as revealed 
in Fig. 5. 

Now we are ready to exhibit that the braiding of edge 
MFs obeys the non-Abelian statistics. As illustrated in 
Fig. 6 with five bricks and three edge MFs, starting from 
the same initial state, the two processes of different orders 
of braiding generate two different final states, which cor- 
responds directly to the inequality UltUtr UtrUlt^ 
namely the non-Abelian feature of braiding. 



setup based on edge MFs than an implementation using 
core MFs simply because the edge MFs are distributed 
over the sample edge. Quantitatively, the excitation gap 
is about O.OlAo for parameters Aq = 0.5to, Vz = O.Sto 
and Or = OMq. For a superconductivity gap Aq ~ ImeV 
[29I . the corresponding temperature of the excita- 
tion gap will be ~ lOOmK, which is not hard for low- 
temperature technologies in these days. 

While not appearing apparently in the criterion for 
chemical potential Vz > \/A^ + /? supporting the topo- 
logical phase, the spin-orbit coupling aR does play an 
important role in determining the magnitude of effective 
p-wave pairing and thus the excitation gap. Besides In As 
[l5| , the layered polar semiconductor BiTel with giant 
spin-orbit coupling can be a promising material for our 
device [37]. The Zeeman splitting Vz of order of several 
meV can be realized in the thin film of strong ferromag- 
netic insulator [38] according to a recent work [29]. 

The present setup has good scalability, and an array 
of the units in Fig. 1(a) supports a ground state wave 
function consisted of a linear combination of the many- 
body Majorana edge states. The ground state can be 
rotated in the degenerate subspace by exchanging pair- 
wisely the Majorana edge states, which is given by a mul- 
tidimensional unitary matrix representation of the 2D 
braid group, and governed by the non-Abelian statistics. 

To conclude, we have demonstrated that the edge MFs 
induced by vortices in topological superconductor can 
be manipulated easily by application of gate voltages at 
point-like constriction junctions. Adiabatic braidings of 
edge MFs have been simulated by TDBdG method, and 
the non-Ableian statistics is proved. The present pro- 
posal therefore provides a promising way for implement- 
ing quantum topological computation based on zero- 
energy MFs, and has the advantage of easy operation 
and simple device structure. 

Methods 

In this work, we use TDBdG method to simulate the 
adiabatic braiding of edge MFs upon switching on 
and off the constriction junctions. The main task 
of the simulation is to solve the TDBdG equation 
(6). The formal solution of the equation is |^(t)) = 
ex.p{—ih J H{t)dt)\'^{0)). To solve it numerically, we first 
divide the total simulation time into small steps of At, in 
which the Hamiltonian can be considered as time inde- 
pendent, and thus |^(t + At)) exp{-ihH{t)At)\'^{t)). 
The exponential can be expanded by the Chebyshev poly- 
nomials Tn{H) 



-ihH{t)At 



(11) 



Discussions 

In an implementation of topological quantum computa- 
tion based on the zero-energy MF modes, the working 
temperature is limited by the energy gap to the lowest 
excitation. The excitation gap is smaller in the present 



Because the coefficients Cm decrease exponentially [33|, 
34], only a finite number (Mmax) of terms in the above ex- 
pansion are needed for sufficient accuracy. Using the re- 
cursive relations of the Chebyshev polynomials Tm{H) = 
2HTm-i{H) — Tm-2{H), the summation in Eq.(ll) can 
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be performed by M^ax times matrix-multiplication. For 
a sparse matrix H in our system, the computation time 

is of o(A^) m. 

The time-dependent part of the Hamiltonian is the 
hopping integrals at the sites on constriction junctions 
between bricks. In the present study, we use a parameter 
< A(t) < 1 to scale the hopping integrals between zero 
to to. In order to simulate adiabatic processes, the func- 
tion X{t) should vary slowly with t, and thus the total 
evolution time T is taken sufficiently long as compared 
with the inverse of the lowest excitation energy. Further- 



more, we use A(t) = (t/T)^ rather than a linear function 
of t, which suppresses efficiently excitations from zero- 
energy MFs to states with finite energy. 
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